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Abstract. The effect of the magnetic field on the en- 
ergy loss rate in the direct Urea reactions is studied. 
The general expression for the neutrino emissivity at ar- 
bitrary magnetic field B is derived. The main emphasis is 
laid on a case, in which the field is not superstrong, and 
charged reacting particles (e and p) populate many Lan- 
dau levels. The magnetic field keeps the process operative 
if Ak/kpn < N F ^ 3 {Np p is the number of the Landau 
levels populated by protons and Afc = kp n — kp p — kp e ), 
that is beyond the well-known switch-on limit in the ab- 
sence of the field, Afc < 0. Cooling of magnetized neutron 
stars with strong neutron superfluid in the outer cores 
and nonsuperfluid inner cores is simulated. The magnetic 
field near the stellar center speeds up the cooling if the 
stellar mass M is slightly less than the minimum mass 
M c , at which the direct Urea reaction becomes allowed 
for B = 0. If B = 3 ■ 10 16 G, the affected mass range is 
M c ~ M < 0.1 M c , while for B = 3 • 10 15 G the range 
is M c — M < 0.015 M c . This may influence a theoretical 
interpretation of the observed thermal radiation as illus- 
trated for the Geminga pulsar. The case of superstrong 
magnetic fields (B > 10 18 G), such that e and p populate 
only the lowest Landau levels is briefly outlined. 



1. Introduction 

The presence or absence of the direct Urea process (n — ► 
e + p + v e , e + p —>■ n + v e ) in the core of a neutron star 
is the most important issue of the stellar cooling. If oper- 
ative, it dominates the cooling at the neutrino stage (age 
t < 10 5 - 10 6 yr) being several orders of magnitude more 
efficient than any other neutrino emission process (e.g., 
Pethick, 1992). However, the direct Urea can occur under 
stringent conditions: one requires quite high fraction of 
protons [kpn < kp p + kp e , where kp a = (37r 2 n Q ) 1 / 3 is a 
Fermi momentum, and n a is a number density of particles 
of species a] in order to conserve momentum in the reac- 
tion. Nevertheless, some equations of state (EOSs) allow 
for that (Lattimer et al., 1991). On the other hand, most 



realistic EOS of dense matter (Wiringa et al., 1988) pre- 
dicts too small fraction of protons, which decreases with 
growing density; consequently, the direct Urea is forbidden 
in the entire neutron star core. 

In this paper we study the possibility for the direct 
Urea to be open in the presence of a magnetic field B, 
if the proton fraction is too low to open the process at 
B = 0. The beta-decay and related reactions in strong 
magnetic fields have been studied since late 1960's (e.g., 
Canuto & Chiu 1971, Dorofeev et al. 1985, Lai & Shapiro 
1991, and references therein). However, these results have 
been obtained under various simplified assumptions (con- 
stant matrix elements, non-degenerate nucleons, etc.) and 
do not give the emissivity of the direct Urea reaction in 
the neutron star cores. Several works on the subject have 
appeared most recently. Leinson & Perez (1997) consid- 
ered the case of superstrong fields (B > 10 18 G), in which 
electrons and protons occupy only the lowest Landau lev- 
els. They found that such fields relaxed the requirement of 
high proton fraction, and the direct Urea was always per- 
mitted leading to a rapid cooling of a neutron star. The 
case of superstrong fields was examined also by Bandy- 
opadhyay et al. (1998). These authors found, by contrast, 
that the condition ll kp n < kp p + kp e " (quotation implies 
that one should be careful with definition of the Fermi mo- 
mentum in superstrong fields) still determined the direct 
Urea threshold, and that the fields enhanced the neutrino 
emissivity by 1 - 2 orders of magnitude in the permitted 
regime compared to the standard value (Lattimer et al., 
1991) 



Q«= 4.10 27 K/n ) 1 / 3 T 9 6 



erg cm s 

m n m p 



(1) 



m* and m* are nucleon effective masses in dense 
and m„ are their bare masses. 



Here Tg — T/10 K, n e is the electron density, hq — 0.16 
fm~ 3 

matter, Tn n oj±±v^ ii^ p 

The present paper is organized as follows. 

In Sect. 2 we obtain a general expression for the neu- 
trino emissivity Q v of the direct Urea process. 

In Sect. 3 we concentrate on a realistic case, in which 
the magnetic field is not extremely high (although still 
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high: B < 3 • 10 16 ), and charged particles populate many 
Landau levels. We show that the field keeps the direct Urea 
open slightly outside the region where /cf« < k Fp + k Fe 
(the standard, B = 0, definition of k F applies). 

In Sect. 4 we illustrate this result by a series of cool- 
ing simulations of magnetized neutron stars with no su- 
perfluid in the inner cores and neutron superfluid in the 
outer cores. 

Finally, in Sect. 5 we treat briefly the case of super 
strong magnetic fields. We show that the results by Lcin- 
son & Perez (1997) are basically correct, although not very 
accurate in details, while those reported by Bandyopad- 
hyay et al. (1998) are inaccurate. 

2. Quantum formalism 

In this section we obtain a general formula for the neu- 
trino energy emission rate in the direct Urea reactions 
valid at any magnetic field B <C I0 20 G (at higher fields 
protons become relativistic) . The calculation is done in 
the Born approximation within the standard quantum 
mechanical framework with weak interactions described 
by the Weinberg-Salam-Glashow theory. We adopt the 
conventional assumption that reacting electrons are rel- 
ativistic, while protons and neutrons are nonrclativistic. 
All these particles are strongly degenerate. The rate of 
transition from an initial state |i) to a final state |/) is 
Q fi = 2iTh~ 1 V- 1 \(f\n\i)\ 2 S{£f - £i), where V is the nor- 
malization volume and Ti describes weak interaction in the 
second-quantization formalism. For the neutron decay, we 
have (c = % = 1) 



H 



G_ 

V2 



(V) 



— QA&tiiGi)^>n V'eT^l + 7 5 )^f(2) 



In this case, G = G F cos9 c , G F = 1.436 • 1CT 49 erg cm 3 
is the Fermi weak coupling constant, 6c ~ 13° is the 
Cabibbo angle, gA = 1.261 is the axial-vector coupling 
constant, Ui is a Pauli matrix, and 7 M is a Dirac matrix. 
Finally, ip a is a field operator in the coordinate representa- 
tion, i.e., the expression of the form ^2 q a q ip q (r), where q 
denotes full set of quantum numbers, d q is an annihilation 
operator, and ip q (Y) is an eigenstate. 

To evaluate the total neutrino energy loss rate (emis- 
sivity) we need to sum Q ^ times the energy of the newly 
born antineutrino over all initial and final states. First of 
all, we can sum the matrix element | (/ 1 Ti \ i) \ 2 over the elec- 
tron spin states. The energy-conserving delta-function is 
not affected by this summation, since all, except the low- 
est, electron states are spin-degenerate. Choosing the Lan- 
dau gauge of the vector potential A = (—By, 0, 0) and 
performing a tedious calculation, we get 
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In these equations, k a i is a cartesian component of a par- 
ticle momentum, e a is a particle energy, s p = ±1 and 
s n = ±1 are, respectively, the doubled proton and neu- 
tron spin projections onto the magnetic field direction, 
k e j_ = \/26n, n is the electron Landau level number, 
b = \e\B, kj, is the antineutrino wave vector, and q = 
(k n x — k vx , k ny — k uy ,0). L x and L z are the normalization 
lengths. Finally, F = F n ^ n [u) and F' — F n i^ n _i(u) are the 
Laguerre functions (e.g., Kaminker & Yakovlev, 1981), n' 
is the proton Landau level number, and u — q 2 /(2b). If 
any index (n or n') is negative, F n ^ n i(u) = 0. 

The next step consists in integrating over the x compo- 
nents of proton and electron momenta, which specify the 
y coordinates of the Larmor guiding centers of these par- 
ticles. This operation gives the factor L 2 Lyb/(2n) 2 and 
removes the second delta-function in Eq. (3). Thus, we 
may write a general formula for the neutrino emissivity 
(including the inverse reaction which doubles the emission 
rate) as 
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dk n dk^ dkp Z dk ez 



x fn (1 - fp) (1 - fe) 5(s n -e e -e p - e v ) e u 
x 3{k n z ~ kez — k pz — k vz ) M. 



(5) 



In this case, f a = (1 + cxp \(e a — ^ Q )/T]) _1 is a Fermi- 
Dirac distribution, and the particle energies are given by 
the familiar expressions: 
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(6) 



with the proton and neutron gyromagnetic factors g p = 
2.79 and g n — —1.91. In principle, the factors gA, g P , g n 
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can be renormalized in dense matter which we ignore, for 
simplicity. 

Since the electron and proton distributions are inde- 
pendent of signs of k ez and k pz we can simplify the ex- 
pression for M by omitting the terms which would anyway 
yield zero after the integration: 



M 1 r /, 2^ 

Q2 = 2 SpS " + 9A > 



1 



h h '. 
• \ p/2 



Thus the situation may be referred to as quasiclassical, 
and corresponding techniques apply. If the main contribu- 
tion comes from large n and n', the difference between F 2 
and F' 2 can be neglected. Moreover, we can neglect the 
contributions of magnetic momenta of particles to their 
energies. Thus, we can pull all the functions of energy out 
of the sum over s 
explicitly: 



and s p , and evaluate the latter sum 



+ 1 



J^M = 2G 2 (l 



■39 a) 
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Inserting this into Eq. (5), and integrating over orienta- 
tions of neutrino momentum we get 
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dkp Z dk ez E v 
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x S(£n ~ £e ~ £p ~ £v) &{k nz — k ez ~ kpz), 

where 



Keeping k vz in the z component of the momentum con- 
serving delta-function in Eq. (5) would lead to a subtle 
thermal effect: it would mollify the resulting functions on 
a temperature scale. We will not pursue the accurate de- 
scription of the effect here, both because the calculation 
would be quite complex, and because the temperature 
scale is assumed to be small. Therefore, we will neglect 
the neutrino momentum in the delta-function and, for the 
same reason, omit it from the definition of the vector q. 
Then, using the isotropy of the neutron distribution, we 
can further simplify the expression for M: 



k 2 

^pz 



#1 



2m; 
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(10) 



(11) 



k p ±_ = \Jlbn' . 
Eq. (6). 

If the magnetic field is not too large the transverse 
electron and proton momenta, k 2 ± and k 2 ± , are sampled 
over a dense grid of values, corresponding to integer in- 
dices n and n' . Thus, the double sum in Eq. (10) tends to 
the double integral over k 2 ± and k 2 ± , 



while the electron energy is still given by 
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G 2 



= 2g\ (8 SpA 5 s ^ 1 F 2 + S Sp ^ 1 S SnA F> 2 ) 26 ]T F„ 2 , » . . . - J 
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(12) 



+ £ Wl+fli) (F' 2 +F 2 ) 
+ S SpSn g A sp (F' 2 - F 2 ) , 



where T represents the small-6 asymptote of the function 



(8) 



where the functions F and F' depend now on u = 
(k 2 nx + k 2 ny )/(2b)^k 2 J(2b). 

Notice that the results of this and subsequent sections 
are equally valid for direct Urea processes involving hy- 
perons. The results for hyperons are easily obtained by 
changing the values of reaction constants (gA, etc.) as de- 
scribed, for instance, by Prakash et al. (1992). 

3. Quasiclassical case 

(a) General treatment and limit B — > 0. First of all con- 
sider the most realistic case of not too high magnetic fields, 
in which electrons and protons populate many Landau 
levels. In this case the transverse wavelengths of electrons 
and protons are much smaller than their Larmor radii. 



26 



F 



Np.Ne 



with N„ 



(«). 



k 2 v J(2b) and N e = k 2 ± /(2b). 



(13) 



After replacing the double sum by the double inte- 
gral Eq. (10) can be considerably simplified. Note, that 
dk z dfc 2 = dk/n = 2m* dek sin#d#, where 9 is a pitch- 
angle. Then the energy integral J de„ de p de e de u is taken 
by assuming, that the temperature scale is small and pro- 
vides the sharpest variations of the integrand. If so, we 
can set k = k F = (S^n) 1 ^ 3 , e = Sp etc. in all the other 
functions. In principle, this assumption constraints the va- 
lidity of the quasiclassical approach. We will come back to 
this point at the end of this section. Finally, we integrate 
over the azimuthal angle of the neutron momentum, and 
over its polar angle to eliminate the momentum conserving 
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delta function, and obtain: 



Qv — Q v x Rg j 

45 7 7 rG 2 (l + 3 fl 2 ) 
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X 0(kFn — \kFp COS p + kp e COS e 



(14) 



where Q® is the field-free emissivity (1), and the factor 
R B C describes the effect of the magnetic field. k p , e ± = 
kFp,e sin Op^e, and k nJ _ is now given by k Fn — (kpp cos P + 
kFeCOsOe) 2 ; 9{x) = 1 for x > 0, 0(x) = for x < 0. 

The asymptotic behaviour of (13) depends on the re- 
lation between the argument of the function F and its 
indices. In the small-6 case one can distinguish three do- 



mains of these parameters: (I) k n j_ < \k p 



I; (ii) 



\k p ± -k e ±\ < k nA _ < (k pA _ + k eA _); and (III) (k pl _ + k eA _) < 
k n ±. In domains (I) and (III) the asymptotes in question 
decay exponentially when k n ± departs from the domain 
boundaries (which are the turning points of corresponding 
quasiclassical equation, e.g., Kaminker & Yakovlev, 1981). 
In both cases the exponents are inversely proportional to 
b, and both asymptotes tend to zero as b — > 0, although 
nonuniformly in the vicinities of the turning points. In do- 
main (II), the asymptote of the expression (13) oscillates 
according to 

^«^cos 2 $, p = \J ANpN e - (Np + N e - uf , (15) 

with a prefactor that is actually independent of b. The 
cosine phase is 



$ = -(1 + N p ) arg(u + N p - N e , -ip) 

+ N e aig(u + N e - N p , ip) - , 

Ja = P(N P -N e ) + i [{Np - N e ) 2 - u(N p + JV e )] 
2u^/N~N~ e 



(16) 



where arg(z) is an argument of complex number z which 
lies in the range [— 7r,7r]. For integer indices, this for- 
mula coincides with Eq. (30) in Kaminker & Yakovlev 
(1981), and it provides an accurate extension of T to 
non-integer indices. One can easily show that a natural 
continuation of Fj ip Ne {u)/{2b) to this case is given by 
the analytical function W^ m (u)/[2buT(N e + 1)T(N P + 1)], 
where Wk m (u) is the Whittaker function, u = k 2 j_/(2b), 
k = (l+i\r p +iV e )/2, and m = (N e -N p )/2. Equations (15) 
and (16) are derived using accurate small-6 asymptotes of 
Wkm{u), which, in turn, could be derived from the inte- 
gral representation of this function by the saddle-point 
method. 

This information is sufficient to check, if R B C repro- 
duces the well-known step-function Q(kp P + kp e — km) 
in the b — > limit. If kFn > kp P + kp e , we are always 



in domain (III) (since the z-momentum is conserved) , and 
R B C — > 0. If, on the contrary, kp n < kp p + kFe the integra- 
tion in (14) covers all three domains, and the boundary of 
domain (II) corresponds to vanishing square root in Eq. 
(15). The contribution from domains (I) and (III) is again 
zero, while the integration over domain (II) yields exactly 
1. To verify this, one can put the rapidly oscillating cos 2 $ 
equal to 1/2 in Eq. (15), and change the integration vari- 
ables: s = cos Op + cos0 e , t = cos Op — cos0 e . In all cases, 
theta-function in Eq. (14) plays no role. 

(b) Effect of the magnetic field in the forbidden do- 
main (kF n > kp p + kp e ). In a finite magnetic field, the 
border between the open and closed direct Urea regimes 
is expected to be smeared out over some scale depending 
on the field strength. This can be important for neutron 
star cooling, since the direct Urea process can remain the 
dominant energy loss mechanism even if it is suppressed 
exponentially by several orders of magnitude. 

To investigate this possibility, one has to calculate 
the emissivity Q v from Eq. (10). We have used two ap- 
proaches, one of which is essentially quasiclassical and 
good for a rapid computation, while the other is of quan- 
tum nature — more precise and time-consuming. 

At the moment it is convenient to introduce two pa- 
rameters, x and y, which characterize the reaction kine- 
matics with respect to the magnetic field strength: 



k 2 



(k Fp + k Fe ) 2 



kl„N 



-2/3 



N 



Fp 



2/3 



(17) 



Fp ly Fp 



where Np p = kp p / (2b) is the number of the Landau lev- 
els populated by protons. The most interesting situation, 
from practical point of view, occurs if x is positive and not 
very large (say x < 10), since whenever kFn significantly 
exceeds kFp + kpe, the reaction is suppressed too strongly. 
The first approach is based on Eq. (14). The main con- 
tribution to this integral comes from the vicinity of the 
line Op — e , corresponding to the least distance from k n ± 
to domain (II). Since this distance should not be large, 
we may use in (14) the small-& asymptotic form of (13) 
in the neighbourhood of the right turning point, i.e., for 
k n ± w k p ± + k e ±. It is given, for instance, by Kaminker & 
Yakovlev (1981) and reads 



^Np,Ne( u ) 



(2b)^(k p± k e± )-^ A;2 



(k p± + k e± )V3 
(smOpSinOe)- 1 / 3 
y (sin p + sin^e) 



2/3 



Ai^) 
Ai 2 (£), 



(18) 
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/<->J.\2 /3 1 ;„ i 7, \4/3 
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- 2y (1 - cos {Op 6 e ))] \ p ■ e > 
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where Ai(£) is the Airy function, defined as in Abramowitz 
& Stegun (1970), and we have assumed kFp — k Fei which 
implies absence of muons and hyperons in neutron-star 
matter. At negative arguments, which correspond to do- 
main (II), the Airy function oscillates, while for positive 
arguments, in domain (III), it decreases exponentially ap- 
proaching the asymptote 

Ai2 «>-i4 cxp (-H' (>>L (19) 

Inserting the latter equation into Eq. (14), we find, that 
at relatively large x 

Rb~xI ^ttt- "4? expf-^V (20) 

At very large x, however, this asymptote should not be 
taken too literally because Eq. (18) ceases to be valid far 
from the turning point. 

On the other hand, under quasiclassical assumptions, 
we are always interested in the case of y ^> 1. If so, we can 
further simplify Eq. (14) at any x using (18), and noting 
that all the functions can be expanded around the line 
9 p = 6 e . Then we obtain 

/OO />7T 
ds / d6sm 2 / 3 9 Ai 2 (0, (21) 
-oo JO 

_ x + s 2 

4 ~ 2 4 /3 sin 2/3 0- 

At x = this integral is taken analytically and gives R B C = 
1/3. By inserting Eq. (19), one easily verifies, that Eq. (21) 
reproduces also the asymptote (20) for x -C 12y. 

Finally, we have calculated the factor i?^, using both 
prescriptions (14) + (18), and (21) at N Fp = y 3 / 2 = 
100, 1000, 10000 and x from to 20. In both these cases 
we have obtained identical results. This indicates that for 
such combinations of x and y, the y-dependence of i?^ c is 
insignificant. These results are plotted in Fig. 1 by open 
circles. The short-dashed curve represents the asymptote 
(20). 

In the quasiclassical approach we have made two ap- 
proximations: firstly, we have replaced the sum by the in- 
tegral at finite 6, and, secondly, we substituted the asymp- 
tote in the form of the Airy function for the function F. 
To assess the quality of both assumptions we have used a 
quantum approach based directly on Eq. (10). 

In this approach, if one transforms the integration over 
k z to that over e in a straightforward manner, the in- 
tegrand becomes singular (the denominators of the form 
^Je — e n appear). If T — > the quantity Q„/T 6 diverges 
at integer Np p . For nonzero T this quantity remains con- 
vergent but oscillates as a function of B and/or density. 
These oscillations are quite familiar and appear in many 
studies (magnetization, electrical and thermal conductivi- 
ties, etc., Landau and Lifshitz, 1986). They are associated 
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Fig. 1. Logarithm of Rb as a function of x in the forbidden 
domain Ak = tf„ — kF V — kFe > for various values of Nf p - 
Open circles correspond to the quasiclassical approach (-R^ c ), 
insensitive to the value of Nf p , while solid circles show the 
quantum numerical results. The short-dash line is the asymp- 
tote (20). Other curves are calculated from the fitting formula 
(22). 



with population of the Landau levels by electrons and pro- 
tons due to variation of plasma parameters. 

If T is larger than the energy spacing between the Lan- 
dau levels the oscillations are washed out, and a smooth 
curve emerges. This regime requires very accurate inte- 
gration over particle momenta in order to include ther- 
mal effect. We were able to perform it only for rather 
small Np P < 20 and do not report these results here. In a 
more important case of lower temperatures the summation 
over Landau levels and energy integration arc indepen- 
dent. The actual neutrino emissivity does oscillate but the 
quantity of practical significance is the emissivity averaged 
over the oscillations (a smooth curve again). We call this 
approach quantum since it involves the summation over 
Landau levels explicitly. In this way we have calculated 
the emissivity and smoothed the oscillations artificially by 
two different methods. We have checked that both meth- 
ods yield nearly identical smoothed Q v . Simple consider- 
ation (see the end of this section) shows that this (non- 
thermal quantum) approach is valid for T < 3 W 'Npp^B, 

or, equivalently, Np p 3 < Hp/T, where is the proton 
gyrofrcquency and \x p is the proton chemical potential. 

The results of these calculations are presented in Fig. 
1 for N Fp = 100, 200, and 400 by solid circles. It is seen 
that with increasing Np P the quantum factor Rb tends to 
the quasiclassical one i?^. 
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We have also found the fit expression that describes 
accurately (Fig. 1) the quantum calculations for < x < 
20 and Nf p > 100 and reproduces the quasiclassical curve 
in the limit Np p — > oo: 



Rr — 



(3x + 6.800) R 



exp 



(x c + 6.800)(3 + xVl2) 

R = Pi + P2C PsX + P^x, 

3.801 2.242 



(22) 



Pi = 1 + 



A7-0.5280 - 
Iy Fp 



Pi = 1 + 



/yO.3498 
JV Fp 



-Pi, 



0.07221 + 



1.751 

atO.5001 ' 
Iy Fp 



Pa = 



1.551 

7\r0.8128 ' 



= xVx + 0.4176 - 0.04035: 



(c) Effect of the magnetic field in the permitted do- 
main (kp n < kp p + kp e ). Since Rb < 1 at s = 0, one 
may expect that the magnetic field has a non-trivial ef- 
fect on the neutrino emissivity in the permitted domain. 
This appears to be true. Applying the quasiclassical ap- 
proach in the form of Eq. (21) at negative x, we obtain for 
R B C an oscillating curve, shown in Fig. 2 by the open cir- 
cles. Using the more accurate quantum approach we get a 
series of the curves for Nf p = 100, 200, and 400, the curve 
with highest Nf p being again the closest to the quasiclas- 
sical result. These oscillations are of quasiclassical nature 
and have nothing in common with the quantum oscilla- 
tions discussed above. From the practical point of view, 
they are not very important, as they hardly have any no- 
ticeable effect on the neutron star cooling. Note, that the 
results for —20 < x < and Np p = oo are accurately 
fitted by the expression (solid curve in Fig. 2) 



Rb 



= 1 - 



COS if 



(0.5816+ IxQ 1 - 192 ' 
1.211 + 0.48231x1 + 0.8453|x| 2 - 533 



(23) 



1.4381a: 



1.209 



To summarize, we remind that the direct Urea reac- 
tion in the B = case is operative if the momentum excess 
Ak = kp n — kpp — kp e < 0. In the presence of the mag- 
netic field, the condition becomes less stringent, and the 
reaction becomes quite efficient as long as x < 10, i.e., 



< N; 



-2/3 



If, for instance, B = 10 16 G, and the 



Ah/ k Fn -i. fp 
density of matter is near the direct Urea threshold, one 
typically has Np p ~ 300 and Ak/kp p < 1/25. 

Notice that the field-free direct Urea process can also 
be allowed beyond the domain Afc < due to the thermal 
smearing of the Fermi surface. If B = and Afc > the 
reaction rate can be written as Q v = x Rp, where 
Rt may be expected to be cx exp(— vp p Ak/T). Thus the 
thermal effect extends the reaction to the domain where 
Ak/kp p < T/n p . The smearing is clearly determined by 
the proton degeneracy parameter, H P /T, which is typically 
about 300 for T ~ 10 9 K. The magnetic field effect is more 
important than the thermal effect if N F ^ < n P /T that can 
often be the case in the inner cores of neutron stars. 
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Fig. 2. Factor Rb in the permitted domain Afc < as a 
function of x for various values of Nf p - Open circles correspond 
to the quasiclassical approach (21), solid line is calculated from 
the fitting formula (23), and the other curves represent the 
quantum numerical results. 



4. Cooling of magnetized neutron stars 

In this section we illustrate the above results by cooling 
simulations. Let us consider a set of neutron star models 
with fixed equation of state but varying central density 
and magnetic field strength. Specifically, we take the EOS 
by Prakash et al. (1988) with the compression modulus 
Kq = 180 MeV and the symmetry energy in the form 
suggested by Page & Applegate (1992). It is assumed that 
matter in the stellar core consist of neutrons, protons and 
electrons (no muons and hyperons). The EOS predicts the 
monotonous increase of the proton fraction with increasing 
mass density. Therefore, if the central density p c is higher 
than the certain threshold density p cr n (corresponding to 
proton fraction x p = n p /rib = 1/9), or, cquivalcntly, the 
stellar mass M is higher than the certain threshold mass 
M c the direct Urea process becomes allowed in a central 
kernel of the inner stellar core. For the chosen EOS, the 
threshold parameters are p cr i t = 12.976 • 10 14 g cm~ 3 and 
M c = 1.442M . 

As pointed out by Page and Applegate (1992) the cool- 
ing history of a neutron star in the field-free regime is ex- 
tremely sensitive to the stellar mass: if the mass exceeds 
M c the star cools rapidly via the direct Urea process, while 
the cooling of the low-mass star is mainly due to the mod- 
ified Urea process, and, therefore, is strongly delayed. The 
effect of the magnetic field would be to speed up the cool- 
ing of the star with a mass below M c , because the strong 
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field opens the direct Urea process even if the B = con- 
dition fcpn < kF P + kFe (or x p > 1/9) is not reached. 

The magnetic field in the neutron star core may evolve 
on cooling time-scales. This may happen if the electric 
currents supporting the field are located in those regions 
of the core, where protons as well as neutrons are non- 
superfluid. If so, the electric currents transverse to the 
field may suffer enhanced ohmic decay due to magnetiza- 
tion of charged particles (e.g., Haensel et al., 1990). The 
consequences of the decay would be twofold. Firstly, if the 
strong field occupies a sufficiently large volume of the core, 
the decay produces an additional heating, which would 
delay the stellar cooling. Secondly, the field decay would 
reduce the direct Urea losses in the forbidden domain 
decreasing the factor Rb- If, however, the neutrons are 
strongly superfiuid, the enhanced decay is absent (Haensel 
et al., 1990; 0stgaard & Yakovlev, 1992), and the ohmic 
decay time of the internal magnetic field is typically larger 
than the Universe age (Baym et al., 1969). 

On the other hand, microscopic calculations of super- 
fluid neutron gaps in the neutron star core suggest that 
the critical temperatures are rather high at not too large 
densities (say, below 10 15 g cm -3 ), but decrease at higher 
densities (see, e.g., Takatsuka & Tamagaki, 1997, and ref- 
erences therein). Thus the electric currents could persist 
in the outer stellar core, where the neutron superfiuid is 
available, and the enhanced field-decay mechanism does 
not work, while the direct Urea reactions operate in the 
inner core and are not subject to superfiuid reduction. 
This latter scenario we adopt. We assume the presence of 
the magnetic field B in the stellar kernel, where the direct 
Urea process can be allowed, and assume nonsuperfluid 
protons in the neutron star core. Thus, the entire core is 
nonsuperconducting, and the magnetic field does not vary 
over the cooling time (age t < 10 7 yr) being frozen into 
the outer core. 

Let us analyse the cooling stage at which the neutron 
star interior is isothermal (t > 10 2 - 10 3 yr). The surface 
temperature T s (seen by a distant observer, i.e., the gravi- 
tational redshift included) is then determined by the heat 
transport through the neutron star crust and is related 
uniquely to the internal temperature. We will use the cool- 
ing code described, for instance, by Levenfish & Yakovlev 
(1996), and Yakovlev et al. (1998). The effects of General 
Relativity are included explicitly. The neutrino luminosity 
is produced by the standard neutrino reactions in the en- 
tire stellar core complemented by the direct Urea process 
in the inner core. The effects of the neutron superfluid- 
ity on the neutrino reactions and neutron heat capacity 
in the outer core are taken into account as prescribed by 
Levenfish & Yakovlev (1996). In addition, we include the 
neutrino emission due to triplet-state Cooper pairing of 
neutrons (Yakovlev et al., 1998). The dependence of the 
surface temperature on the internal stellar temperature 
is taken from Potckhin et al. (1997) assuming no enve- 
lope of light elements at the stellar surface. To emphasize 



the effect of the internal magnetic fields on the neutron 
star cooling we neglect the presence of the surface mag- 
netic fields and assume that the inner field does not affect 
the relationship between the surface and internal temper- 
atures. 

The density dependence of the neutron critical tem- 
perature (triplet-pairing) is given by the step-function: 
T cn = 10 10 K at p < 7.5 • 10 14 g cm~ 3 , and T cn = at 
p > 7.5 • 10 14 g cm~ 3 The resulting cooling curves, T s (t), 
are insensitive to the initial inner temperature provided 
the latter is sufficiently high ( > 10 9 K). 

The cooling curves are depicted in Figs. 3 and 4. Dash 
line illustrates fast cooling due to direct Urea process al- 
lowed in a large portion of the core at B — 0. While cal- 
culating this curve we have used exact factors Rb in all 
the neutron-star layers where the field-free direct Urea is 
either permitted or forbidden. We have verified that the 
results are insensitive to specific values of Rb in the per- 
mitted domain. In that domain, one can safely use the 
quasiclassical fit (23) or even set Rb = 1. On the other 
hand, even a huge internal field is unimportant in the for- 
bidden domain (the curves for B = and 3 • 10 16 G co- 
incide): new regions of the core, where the direct Urea is 
open by the field, amount for a negligible fraction of the 
total neutrino luminosity. 
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Fig. 3. Logarithm of the surface temperature as seen by a dis- 
tant observer as a function of neutron star age. The dash curve 
is for a star of mass M = 1.595M0, well above the threshold 
mass M c , with magnetic field < B < 3 • 10 16 G. The dotted 
and solid curves are for the 1.439M Q star. Bar shows observa- 
tions of the Geminga pulsar (Meyer et al., 1994). 



The upper dotted curves are calculated for the stars 
with masses 1.439 M (Fig. 3) and 1.320 M© (Fig. 4) at 
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B = 0. They represent the slow cooling via the standard 
neutrino reactions (the direct Urea is forbidden) . The solid 
curves illustrate the effect of the magnetic field for the 
stars of the same masses. If the stellar mass is slightly 
(by 0.2 %) below M c (Fig. 3) the cooling curve starts 
to deviate from the standard one for not too high fields, 
B = 3 • 10 14 G. Stronger fields, 10 15 -3 • 10 15 G produce the 
cooling intermediate between the standard and rapid ones, 
while still higher fields B > 10 16 G open the direct Urea 
in a large fraction of the inner stellar core and initiate a 
nearly fully enhanced cooling. If, however, the mass is by 
about 8% below the threshold one (Fig. 4), only a very 
strong field B — 3 • 10 16 G could keep the direct Urea 
process slightly open to speed up the cooling. 
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Fig. 4. Same as in Fig. 3. The dotted and solid curves are for 
the 1.320M Q star. The dash line is the same as in Fig. 3. 

The results indicate that the magnetic field in the very 
central stellar core can indeed enhance the cooling pro- 
vided the stellar mass is close to the threshold mass M c . 
If B = 3 • 10 16 G, the effect is significant in a mass range 
M — M c < 0.1 M c . For lower fields the range becomes 
smaller. If, for instance, B = 3 • 10 15 G, the mass range 
becomes as narrow as M — M c < 0.015 M c . 

Our results can be used for interpretation of obser- 
vational data. By way of illustration, consider observa- 
tions of the thermal radiation from the Geminga pulsar. 
Meyer et al. (1994) fitted the observed spectrum by the 
set of hydrogen atmosphere models. These fits yield rather 
low non-redshiftcd effective surface temperature T e g = (2 
- 3)-10 5 K. Introducing the appropriate redshift factor 
g = yjl — Rg/R ~ 0.8 (R is the stellar radius and R g is 



1.45 - 




1.35 1 1 

12 3 

B 16 [G] 

Fig. 5. The allowed mass range for the Geminga pulsar as a 
function of the internal magnetic field 

the gravitational radius) one gets redshiftcd surface tem- 
perature lgT s [K] = 5.29 ±0.09 (T s = gT cS ). Adopting 
the dynamical Geminga's age t = 3.4 x 10 5 yr we can 
place the Geminga's error bar in Figs. 3 and 4. Let us 
use our cooling model (with possible strong magnetic field 
B near the stellar center unrelated to the much weaker 
Geminga's surface magnetic field). If B = Geminga is 
found between the lines of standard (M < M c ) and fast 
(M > M c ) cooling. It is clear that tuning the mass slightly 
above M c we can force the cooling curve to cross the er- 
ror bar. However, the mass range corresponding to the 
bar width (Fig. 5) would be tiny (about O.OO3M ), as the 
cooling rate is extremely sensitive to M in the domain just 
above M c . The narrowness of the confidence mass range 
makes it fairly improbable that the Geminga's mass lies 
in this range. Accordingly the suggested interpretation of 
the Geminga's cooling is unlikely. The situation becomes 
strikingly different in the presence of the strong magnetic 
field. The field shifts the confidence range of M (faded 
area in Fig. 5) below M c (cf. Fig. 3), where variation of 
T s with M at a given t is much smoother. This widens 
considerably the confidence mass range. At B = 10 16 G it 
is about 0.02 M , while at B = 3- 10 16 G it is about 0.04 
Mq so that the chances that the Geminga's mass falls into 
this range become much higher. This makes the proposed 
interpretation more plausible. 

5. The case of superstrong magnetic fields 

Using the formalism of Sect. 2 let us outline the main fea- 
tures of the direct Urea reaction produced by electrons and 
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protons occupying the lowest Landau levels (with proton 
spins aligned with the magnetic field). Notice, however, 
that one needs superhigh magnetic field, B > 10 18 G, to 
force all electrons and protons into their ground Landau 
levels. Let n e and n p be the number densities of these par- 
ticles. Contrary to the field-free Fermi momenta valid at 
not too high fields and used in Sects. 3 and 4, the limiting 
momenta along the superhigh magnetic field become field- 
dependent, fcg p = 2-K 2 n e ^ p /b, fc° w fi e . The distribution of 
neutrons can be characterized by two Fermi momenta for 
particles with spins along and against the magnetic field 
kFn,s n , kpn,i < kpn,-i- Our starting point is Eq. (8), 
which reduces to 



M 
F 2 



G 2 F 2 



(1 -g A ) 2 S Sn ,i + 2g\ 6s 



(24) 



00 



exp 



k 2 
2b 



Let us insert it into Eq. (5), neglect k vz in the momentum- 
conserving delta function, and integrate over the neutrino 
orientations and over an azimuthal angle of the neutron 
momentum. Then we convert the integrals over k ez and 
k pz into the integrals over particle energies, take the stan- 
dard energy integral, and perform the integration over the 
neutron pitch-angle. We obtain (the inverse reaction in- 
cluded) 



10080 

b v 

X *°*°(l + 3<fl) J- 

+ g 2 A 0(u_ la )e- u -^] . 



m n m p /i e T 
1 



(1 - g A ) &{ui a ) e~ 



(25) 



- + aklf, and a = ±1 



cor- 



where 2bu Sn . a = fcf 
responds to two different reaction channels in which the 
electron and proton momenta along the z-axis are either 
parallel (a = 1) or antiparallcl (a = —1). The step func- 
tions indicate that the channels are open if u s „, Q > 0. The 
channel a = — 1 is always open in npe dense matter with 
the superstrong magnetic field, while the channel a = 1 
is open only if kFn > k p + h®. The latter condition is op- 
posite to the familiar condition kp n < kp p + kp e in the 
field free case. One has exp(— u Sni<x ) < 1 and Q v oc b 2 , 
but one cannot expect Q v to be essentially larger than 
the field-free emissivity Q° as long as B < 10 19 G. Note 
also, that Eq. (25) describes the contribution of particles 
populating the ground Landau levels to the emissivity not 
only in the limit of superstrong fields. For it to be valid at 
moderate fields one should substitute number densities of 
e and p on the ground levels for n e and n p in the definition 
offc e V_ 

Similar result for the superstrong magnetic fields has 
been obtained recently by Leinson and Perez (1997). Nev- 
ertheless, their expression differs from our in several re- 
spects. Most importantly, the authors got (l + g\) instead 



of (1 — 5a) 2 , which substantially overestimates the partial 
rate in the corresponding channel. Much more different re- 
sult under the same assumptions was obtained by Bandy- 
opadhyay et al. (1998). In our notations, the latter authors 
found the emissivity to be proportional to 0(— iinje"" 11 , 
and obtained the spurious enhancement of the neutrino 
emission and additional acceleration of the cooling in a 
superstrong magnetic field. Their condition which opens 
the reaction channel is opposite to the actual one. 

6. Conclusions 

We have considered the neutrino emission produced by 
the direct Urea process in the cores of neutron stars with 
strong magnetic fields. We have derived the general ex- 
pression for the neutrino emissivity [Eqs. (5) and (7)], and 
analysed it (Sect. 3) in the most important case, in which 
the magnetic field is not too strong, and the reacting elec- 
trons and protons populate many Landau levels. We have 
shown that the magnetic field can switch on the direct 
Urea process under the conditions, in which the field-free 
process is strongly suppressed by momentum conserva- 
tion. We have obtained the analytical fit expressions (22) 
and (23) for this enhanced emissivity. We have performed 
a series of neutron star cooling simulations for the model, 
in which the star has a strong neutron superfluidity in the 
outer core but no superfluidity in the inner core. We have 
demonstrated, that the magnetic fields B < 3 • 10 16 G in 
the stellar center can strongly enhance the neutrino lu- 
minosity and accelerate the cooling of neutron stars with 
masses by about 10 % below the threshold mass, minimum 
mass of a star at which the field-free direct Urea process 
starts to be permitted. Therefore, the cooling can be con- 
trolled by the central stellar magnetic fields. The effect 
does not require superstrong fields B > 10 18 G analysed 
by different authors. 
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